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A quantum theory of feedback of bosonic many-atom systems is formulated. The feedback-induced many- 
atom correlations are treated by use of a parameterized correlation function, for which closed equations of 
motion are derived. Therefrom the dynamics of any additive property of the system, i.e., properties derived from 
the reduced single-atom density operator, can be obtained. An example is given that indicates the correlation 
effects of feedback. 
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I. INTRODUCTION 

Feedback is widely used in engineering to control certain 
properties of a system. Motivated by optical homodyne mea- 
surements, a theory of feedback has been formulated for cav- 
ity fields [|J. Many applications of feedback have been pro- 
posed for single quantum systems, such as, besides single- 
mode cavity fields, also for single atoms or generic har- 
monic oscillators. Moreover, extensions to time-delayed non- 
Markovian feedback have been derived [||] and coherent feed- 
back [0| and feedback using weak measurements [^j have 
been studied. Furthermore, alternative approaches to feedback 
have been presented in the context of quantum-state estima- 
tion and quantum control [p|] . 

Besides laser cooling [n] many advances have been accom- 
plished in the last years in cooling atoms to the ultimate limits. 
Among them are Bose-Einstein condensation of dilute gases 
of alkali species [Q, |i||9| [k]] and atomic hydrogen [|H]|, the ex- 
perimental realization of atom lasers [ |T2[ |l3[ |I4] ] and the im- 
plementation of micro- and nanostructured atom-optical de- 
vices Jl5| , [To] p"7| ]. This experimental progress draws attention 
to the investigation of many-atom systems and their conden- 
sate features from a quantum-optical viewpoint. In view of 
the experimental capabilities up to date, it is quite promis- 
ing to engage in developing methods for controlling complex 
many-atom systems, such as produced nowadays in the labo- 
ratory. 

The use of feedback methods acting on accessible macro- 
scopic observables of the system appears particularly desir- 
able. An experimental implementation of it would not require 
the active intervention of the experimenter, but would act as 
a self-automated basic technological tool, such as it is now in 
engineering. Unfortunately, for the case of feedback acting 
on a multi-mode many-atom system, the existing descriptions 
of feedback are in general not applicable: When trying to ac- 
cess arbitrary properties of the system, these methods become 
intractable, due to the required vast number of degrees of free- 
dom. 

Typically employed mean-field approximations, could be 
used to factorize or truncate higher-order many-atom corre- 
lations, in order to simplify the problem. In the considered 
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case, however, the feedback operation contains the measure- 
ment of a macrostate observable of the many-atom system, 
to which all atoms contribute. The collective access to the 
atoms in such a measurement induces strong correlations be- 
tween them. These correlations represent a main ingredi- 
ent of the feedback and should be taken care of in a non- 
approximate way. The formation of correlations becomes 
clear when considering that an observed measurement out- 
come of a macrostate observable can be produced by a vast 
number of different microstates, that all enter the feedback op- 
eration in a coherent superposition. Such superposition states 
may contain strong correlations between the atoms and one 
may therefore expect that in general many-atom quantum cor- 
relations play a dominant role in the dynamics of the feed- 
backed system. 

Correlations of this type are usually generated by interac- 
tions between atoms, such as for example two-body collisions. 
In the case considered here, such explicit interactions will not 
be taken into account. The many-atom correlations will solely 
be generated by the measurement process. An interaction be- 
tween the atoms can in fact be identified in the measurement. 
The measurement process can be thought of as a probe field 
interacting with all the atoms with a subsequent projection of 
the probe by the final readout of the measurement result. Thus 
we may imagine an effective interaction between the atoms 
being mediated via the probe field. 

In the following an exact method is presented for describ- 
ing the reduced single-atom dynamics under the influence of 
feedback operations acting on macroscopic observables of a 
bosonic multi-mode many-atom system. The feedback loop 
consists of a measurement with predefined measurement un- 
certainty (7 and a subsequent unitary operation that again acts 
on all atoms of the system. The feedback loop is assumed to 
act instantaneously without time delay between measurement 
and operation. As will be shown, the many-atom correlations, 
that are generated by the measurement, can be formally incor- 
porated into the single-atom dynamics by introducing a dy- 
namically evolving additional degree of freedom. Thus the 
free evolution intermitted by feedback operations can be for- 
mulated in terms of a quasi single-atom dynamics, revealing 
modifications to the feedback dynamics of a true single-atom 
system. 

The paper is structured as follows: In Sec. II the macro- 
scopic observables needed for the feedback are introduced. 
Furthermore, the effect of the feedback on the many-atom 
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FIG. 1: Schematic function of the feedback loop. First a measure- 
ment on the system occurs, described by the operator M(A— A), The 
value of the measured outcome A is then used as a parameter for the 
unitary operation U(A). 



quantum state is described and a first indication of the im- 
portance of atom-atom correlations is given. Section |l| pro- 
ceeds by introducing a specially designed correlation function 
that is shown to allow for a closed description of the feed- 
back. Besides the dynamical evolution during the free evo- 
lution between feedback operations, also an interpretation of 
this correlation function is presented. An application of the 
developed method to the case of a condensed bosonic gas is 
given in Sec. and finally a summary and conclusions are 
found in Sec. W\ 



II. FEEDBACK LOOP 

The feedback loop shall be considered as a measurement of 
a macroscopic observable A of the system and the use of the 
so obtained information, i.e., the measured value A, for apply- 
ing a unitary operation, U(A), on the system. Information is 
therefore gained by measurements and used to manipulate the 
system, cf. Fig. ||. For the macroscopic observable we may 
think of any extensive quantity of the system, i.e., a quantity 
that has additive contributions from each atom in the system. 
Thus measurements are not performed on single atoms but on 
large numbers of atoms. 

First let us consider the single-atom observable A and its 
canonical conjugate B with commutator (h= 1) 

[A,B] = i. (1) 

Note that from now on for denoting single-atom Schrodinger 
operators calligraphic symbols are used. We now turn to a 
second-quantized description of the many-atom system by in- 
troducing the bosonic atom-field operator 4>{x) with commu- 
tator 

[$(x),ft(x')] = 6(x-x'). (2) 

Using an expansion into a set of complete and orthonormal 
modes u\(x), 

$(x) ='^2<i>xU\(x), (3) 

A 



the resulting commutator relation for the annihilation and cre- 
ation operators of atoms in mode A, <fi\ and (fr x , respectively, 
reads 

[$x,$] i ]=6 ltX . (4) 

The macroscopic observable A that is used for the mea- 
surement in the feedback loop is a sum of contributions of the 
individual atoms. In the formalism of second quantization A 
is therefore bilinear in the atom-field operators and reads 

//A 

Here (fi\A\X) are the matrix elements of the correspond- 
ing single-atom operator A with the single-atom Schrodinger 
modes |A) being defined by: (x\X) = u\(x). 

For implementing a feedback operation U(A) that is ca- 
pable of changing the measured property A, the canonic ally 
conjugate variable of A is needed. Whereas for a single atom 
this is trivially given by the operator B, cf. Eq. ([l]), for a gen- 
eral many-atom system a conceptual problem arises. Clearly, 
the proper canonical conjugate is given by 

& = ^J2ftM»\B\X), (6) 

/^A 

where again (fi\B\X) are the matrix elements of the corre- 
sponding single-atom operator B. Moreover, A is the operator 
of the total number of atoms in the system, 

A 

With this definition the sought commutator relation is estab- 
lished: 

[A,B]=i, (8) 

as in the case of a single atom, cf. Eq. (jl|). 

However, due to the appearance of the atom-number op- 
erator N in Eq. (^), the operator B is not bidiagonal in the 
atom-field operators. It is an intensive quantity and not an ex- 
tensive one and thus is not a simple sum of contributions of 
the individual atoms. Thus it cannot be experimentally real- 
ized by application of an external field (e.g. laser field) that 
interacts with each individual atom. In other words, the prob- 
lem is, that the total number of atoms is in principle unknown 
without its particular measurement. For the feedback action 
we need to know how many atoms have contributed to the 
previously measured quantity A. The lack of this knowledge 
is a potential source of imperfection of the overall feedback 
action. 

This problem can be overcome in two ways: Either in each 
feedback loop not only A but also the atom number A is mea- 
sured, or the atom number is estimated. We will here concen- 
trate on the latter approach, where the estimated number of 
atoms A c may be derived from a priori knowledge about the 
system, without additional measurements. 



3 



It is worth noting, that neglecting this issue by replacing N 
by its average (N) in Eq. (^), but still assuming the commuta- 
tor relation (|8]), results in the neglect of atom-number fluctua- 
tions. Such an approach is only valid when the atom number 
is exact. Considering the case where only part of the system 
is subject to the feedback, say only atoms being located in 
a predefined spatial region, an exact atom number in this re- 
gion is rather doubtful even when the total number is precisely 
known. 

Such an approximation can be avoided by consistently 
defining the operator B by use of the estimated atom number 

N c , 



B = 



and using the resulting commutator relations 

[A, B] = iN/N Cl [A, N] = [B, N] = 0. 



(9) 



(10) 



The measurement outcome is then a pure addition of all 
single-atom contributions. And the subsequent unitary oper- 
ation will feed back the appropriately estimated part to each 
atom. 

Let us now consider the action of the feedback loop on the 
quantum state of the many-atom system. Let t T be times in- 
finitesimally before and after the time t when the feedback is 
applied. The action of the feedback on the many-atom density 
operator g then reads 

g(t+) = J dA U(A) M(A—A) g(t~) M\A — A) U\A). 

(11) 

It expresses the relation between the many-atom density oper- 
ators g(t~) and g(t + ) before and after the feedback, respec- 
tively. The function M(A) is the resolution-amplitude [ plj ] 
that describes in its operator form the measurement process 
with a measurement uncertainty a determined by the width 
of |M(A)| 2 . For simplicity we may take it to be a Gaussian. 
The integration signifies the classical averaging over all possi- 
ble measurement outcomes. Note that Eq. (pTl) can hardly be 
solved due to the large number of degrees of freedom of the 
many-atom density operator. 

An important class of feedback mechanisms is defined by 
the following action of U(A) on the observable A, 

tft(A) Au(A) = A + f(A) N/N. (12) 
It is realized by choosing for the unitary feedback operation 



U(A) = exp -if (A) B 



(13) 



Here the function f(A) determines the response of the feed- 
back loop to a measured value A. It encodes the desired action 
of the feedback and is a numerical function chosen by the ex- 
perimenter. 

An example may illustrate the feedback operation given in 
Eq. (|l3]): For a system containing exactly N atoms (N e = N) 



and choosing for example f(A) — s(A+Aq), the transforma- 
tion ( |l2| ) corresponds to positive (s > 0) or negative (s < 0) 
feedback with an offset Aq. In the case of negative feedback 
and vanishing offset (s = — 1, Aq = 0) we may derive from 
Eq. < Jl~T| ) relations of the single-atom quantities A and B be- 
fore and after the feedback. That is, the system is measured 
and fed back as a whole, and we regard the effects on the av- 
erage properties of an individual atom in the system. 

The relations for the mean values before and after the feed- 
back read 



(A) t+ = 0, (B) t+ = (BY. 



(14) 



Thus on average the quantity A is perfectly compensated to 
zero without affecting B, as expected for this type of feed- 
back. Conceptual differences from the behavior of a true 
single-atom system are revealed when regarding the variances 
of these properties. The variance of B reflects the expected 
increase of noise due to the back-action of the measurement: 



(AB)' 



(AB) 2 



(15) 



However, the variance of A, i.e., of the single-atom version of 
the measured observable, does not correspond to that expected 
for the feedback acting on a single-atom system: 



(AAY 



t+ 



(AAf 



N 2 



a \ 2 



(16) 



For N — 1 atom in the system we have A^ A and there- 
fore the first two terms on the r.h.s. cancel each other, leaving 
as variance only the measurement uncertainty a 2 . This is the 
result that we expect for the single-atom case. In the general 
N atom case, however, the first two terms do not cancel each 
other. Equation ( fTo] ) can be interpreted as follows: The uncer- 
tainty of the single-atom observable A is reduced by the gain 
of knowledge on the macroscopic observable A. However, 
the uncertainty of the single-atom quantity A is in general not 
identical to the uncertainty of the macroscopic observable A 
divided by the number of atoms. And therefore the first two 
terms do not compensate and in this way increase the noise 
above the measurement uncertainty term. 

The second term in Eq. ( |T6| ) contains atom-atom correla- 
tions of the type ($j^A$L$A')' as can be seen from the def- 
inition (|5|). Therefore, the amount of compensation of noise 
in the single-atom properties due to the feedback depends on 
many-atom quantum correlations. The determination of these 
correlations requires a hierarchy of equations for higher-order 
atom correlations. To incorporate all these correlations in a 
consistent way is therefore required even for obtaining single- 
atom properties. 



m. SOLUTION IN TERMS OF CORRELATION 
FUNCTIONS 

Our goal is to describe expectation values of arbitrary ex- 
tensive properties of the system, such as for example the total 



4 



energy. Any such quantity can be expressed in terms of the 
so-called single-atom density matrix p^x, that is defined by 

This density matrix describes the average properties of single 
atoms in the system and is normalized to the total average 
number of atoms 

J2pxx = (N). (18) 

A 

However, for the purpose of including the necessary higher- 
order atom correlations and for obtaining a closed set of equa- 
tions we study a slightly different quantity. Namely, we con- 
sider the change due to a feedback operation of the correlation 
function, 

D MX (z) = (${f(z)$ li y (19) 

The operator T(z) is yet unknown, but we choose it to satisfy 
the requirement 

lim T(z) = 1. (20) 

In this way the single-atom density matrix results as the value 

for 2 = 0, 

D M (Q) = (0t^>. (21) 

As already mentioned, with the introduction of T(z) higher- 
order many-atom correlations are contained in D li \{z), which 
is the clue to the solution. 

From Eq. ( |TT| ) the correlation D^\(z, t + ) right after a feed- 
back operation can be expressed as 

D^(z,t+) = fdA^(A-A)U^A)4 > lf{z)4> fl 

xU(A)M(A-A)^ _. (22) 



This kernel can easily be calculated since it only contains ma- 
trix elements of the single-atom measurement operator, that is 
obtained from the resolution amplitude M(A) by substituting 
for the argument the single-atom operator A — A, i.e., 

M(A-A) = M(A-B)\ B=A . (27) 
Moreover, it contains the single-atom unitary operation that 



Now the operator function T(z) is chosen in such a way, that 
the r.h.s. of this equation can again be expressed in terms 
of the correlation functions. More precisely, in terms of 
Dfj,\(z,t~) at the time t~ before the feedback operation. It 
can be shown that the specific Ansatz 



f (z) = exp(iz T -Z), (23) 



is appropriate for that purpose. Here the vector parameter is 
z T = (a, /?, 7) [ [22I ] and Z is defined as the vector operator 



Z = B . (24) 
V N/N c ) 



Using this Ansatz, from Eq. ( p2[ ) we obtain a closed set of 
equations that connects the correlations _D m a(z) immediately 
after the feedback operation with those before the feedback: 



ZVM + ) = E j d 3 z'F^ x '(z,z')D^x>(z',t-). (25) 



The function (z, z') represents the action of the feed- 
back mechanism and will be denoted in the following as feed- 
back kernel. It is defined as 



(26) 



reads 

U{A) = exp[-t/(A) B/N c ] . (28) 

The kernel (^6|) can be interpreted in terms of single-atom 
properties as follows: First the single atom's observable A is 
measured with outcome A' ± A" and then this observable is 
shifted by f(A)/N c , as would be expected when equally dis- 
tributing the collective shift f(A) among the estimated num- 



i^; V (a,/?, 7 ,a',/?',7') = S(/3-(3') JL J ' d A J dA' J dA" J dk «S[ 7 - 7 ' + af(A)+f3k] 

x exp[i(a-a')(A-A')+ikA"] (^\U(A) M(A' +A" -A)\ f /) (A'| M\A' - A" -A) U\A)\X). 

I 
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ber of atoms. The connection between A and A'±A", how- 
ever, is determined by the mapping of the additionally intro- 
duced degrees of freedom, z and z'. This latter mechanism is 
the way how all orders of many-atom correlations are incor- 
porated in the mapping and signifies the deviation from the 



J 



true single-atom case. 

Considering the case z = (a — (3 — 7 = 0), from Eqs ( |25j ) 
and (|26|) it becomes immediately clear that for the single-atom 



density matrix 



} no closed mapping is obtained: 



1+ 



= Jda'JdAjdA', 



—iot'{A—A') 



{li\U{A)M{A'-A)\n') 



x{\'\M\A'-A)U\A)\\) D^ x ,(a',t-), 
I 



(29) 



where 



D fl - V (a',r) = v(z',r)| 



= (a',0,0)' 



(30) 



On the r.h.s. of Eq. ( [29] ) higher-order many-atom correlations 
always remain in D^x' («', t~) due to the non-vanishing inte- 
gration variable a' . That is, even when considering the single- 
atom dynamics, the mapping of the full correlations D^\(z) 
has to be employed. Afterwards z can be set to zero to obtain 
the sought single-atom density matrix. 

It is worth mentioning, that an approximate closed mapping 
for the single-atom density matrix in fact exists in the limit 
(7 — > 00. Then the measurement resolution amplitude M. only 
slowly depends on its argument and the A' integration in ( p9[ ) 
leads to non-vanishing contributions only for a' « 0. That 
shows that for very bad measurement resolutions, almost no 
correlations are generated between the atoms. 

In the special case when exactly one atom is present in the 
system, the correlation function will naturally reduce to the 
single-atom density operator D ll \(z) — > ((/>[ fy^). Therefore, it 
does no longer depend on z. In this true single-atom case, the 
a' integration in the mapping (|5|) acts only on the kernel and 
results in the usual relation of the form given in Eq. (11). Now, 



however, for the single-atom density operator of a true single- 
atom system. Thus, the need of the parameter z signifies the 
presence of many-atom correlations in the feedback mapping. 

Between two feedback operations the many-atom system 
will evolve according to the Hamiltonian H. Assuming non- 
interacting atoms it is of the form 



ff = 5>|7i|A) $l$ x . 



(31) 



Without loss of generality the mode functions u\(x) — (x\X) 
are chosen as energy eigenstates of the single-atom Hamilto- 
nian H with energies E\. Given the unitary time-evolution 
operator 



U(t) = exp(-iHt), 



(32) 



we may employ the time-dependent version of the operator 

T(z), 



f (z, t) = C/t (t) f (z) U{t) = exp \i z T - Z(t) 



(33) 



where Z(t) reads 



Z(i) = U\t)ZU{t). 



(34) 



Using these definitions we arrive at the expression for the cor- 
relation (z) propagated according to the Hamiltonian (|3lj 
from time t' to t, 



f(z,t-t')< 



-iiE^-Ex^t-t 1 ) 



(35) 



Now an assumption has to be made with respect to the free 
time evolution of the operators contained in Z. We assume 
that Z(t) can be given as 



Z(t)=V z (t-t')-Z(t'), 



(36) 



where the 3x3 real-valued matrix has the property V z (— t) = 
V^ 1 (i). Note that the free time evolution of many physical 
systems can be cast into this form. Representative examples 
are free or harmonically trapped atoms when A and B are 
center-of-mass position and total momentum or vice versa. In 
such cases the correlations evolve as 

D mA (z, t) = Jd 3 z' W„a(z, z', t-t') D mA (z', t'), (37) 



where in energy representation the kernel is given as 



^ A (z,z',<)^ 3 )[z'-Vj(;)-z] 



,-i(E p -E A )t 



(38) 



which again can be evaluated rather easily. 

Now it becomes also clear, why the additional parameter f3 
has been introduced. Clearly it is not necessary for describing 
the effect of a feedback operation on the single-atom density 
matrix, as can be seen from Eq. (|29|). For the free evolution, 
however, this parameter is needed to provide a complete set 
of dynamically evolving variables. Consider, for example, the 
case of harmonically trapped atoms where the feedback is per- 
formed on the total momentum. Then a is associated with the 
total momentum, whereas j3 is associated with the center-of- 
mass coordinate. The two parameters are then coupled during 
the free evolution, which shows the need for the additional pa- 
rameter (3. Note, that the parameter 7 naturally emerged due 
to the commutator relation of A and B. 
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It has been shown that by introduction of an additional 
(real-valued) vector parameter z a specifically chosen corre- 
lation function can be evolved in time by closed sets of equa- 
tions of motion. Both for the feedback of a macroscopic ob- 
servable and for the free evolution between the feedback oper- 
ations [under the constraint (f3o|)] the appropriate kernel func- 
tions can be given. The vector z could be formally interpreted 
as being an additional degree of freedom of the single atom, so 
that we may view the complete feedback dynamics as a quasi 
single-atom problem. This is justified also by the fact that all 
required kernel functions are obtained from matrix elements 
of single-atom operators. 

The correlation function D^\(z) reveals another interest- 
ing aspect. Clearly, for z = it converges to the single-atom 
density matrix and describes the microscopic properties 
of the average single atom. It has however also a connection 
to macroscopic properties, as will be shown now. If we would 
neglect in the definition ( |l9| ) the field operators tp\ and </> M , we 
could approximately write 



AP/Apo 




1 23456789 10 

TV 



FIG. 2: Variance of the single-atom momentum, scaled in units of 
the ground-state variance Ap , versus atom number TV for the state 
after the feedback. Measurement uncertainties are a /Ap = 1 (filled 
boxes) 1.5 (empty boxes), and 2 (triangles). 



D mA (z) ~ D(z) = (f (z)>. 



(39) gas reads 



Given the definition of T(z) in Eq. ( |23| ) it becomes clear that 
D(z) is the characteristic function (i.e. Fourier transform) of 
the probability of measuring simultaneously the observables 
A, B, and TV/TV e . For a system with precise atom number, the 
latter observable can be omitted and we essentially arrive at 
the Fourier transform of the Wigner function W(A, B) of the 
pair of canonically conjugate macroscopic variables A and B. 

Turning back to the true definition (|l9|) it can now be ob- 
served that the diagonal correlation D\\(z) can also be inter- 
preted as the Fourier transform of a Wigner-type function of 
the macroscopic phase space. However, it is that of a many- 
atom state, where one atom in mode A has been removed. 
The off-diagonal elements (z) (A ^ /i) represent the cor- 
responding terms including single-atom coherences. Thus, 
the chosen correlation function covers both microscopic and 
macroscopic properties. Moreover, it accomplishes this in a 
specific combination, that allows for the solution of the feed- 
back dynamics. 



IV. APPLICATION TO AN IDEAL BOSE GAS 

For illustrating the effects of atom correlations that can be 
treated with the presented method, let us consider feedback of 
a gas of TV trapped bosonic atoms. We assume all the atoms 
to be initially in the ground state of the harmonic trap poten- 
tial, e.g., a simple model for a non-interacting Bose-Einstein 
condensate. The feedback shall consist in measuring the to- 
tal momentum P of all atoms and in subsequently shifting 
this observable to zero by use of the observed value and the 
center-of-mass operator Q [A— *P, B^Q, f(P) = —P, and 
N C = N]. 

The initial correlation function of the degenerate bosonic 



D M (z) =JV(/z|0}(0|A)e X p[i 7 (l-l/iV)] 

(Ag ) 



x exp 



TV-1 



2„2 



(A Po ) 2 a 



TV 2 



(40) 



Here Aq and Ap are the position and momentum uncertain- 
ties, respectively, of the single-atom ground state |0) in the 
trapping potential. This example shows quite well the above 
mentioned combination of microscopic and macroscopic de- 
scriptions. In this special case the correlation ( pfj| ) is a prod- 
uct of the single atom density matrix = TV(/i|0) (0|A), a 
phase factor determined by the precise number of atoms, and 
a Gaussian distribution. The latter is the Fourier transform of 
the joint probability distribution for simultaneously observing 
the macroscopic values Q and P of TV — 1 atoms. That is, it is 
the Fourier transform of the Wigner function W(Q, P) of the 
center-of-mass degree of freedom of the macroscopic atomic 
gas with one atom being removed. 

Though the correlation ( |40| ) could be propagated in time 
as described previously, we consider here only the result of a 
single feedback operation. In this way, the variances of single- 
atom quantities, such as position and momentum, can be stud- 
ied in dependence on the atom number TV and the measure- 
ment uncertainty a. Figure || shows the variance of the single- 
atom momentum in dependence on the atom number for three 
different measurement uncertainties a. One observes that for 
larger numbers of atoms, the momentum variance converges 
to the initial ground-state variance Ap , which establishes the 
regime of weak measurements where the single atom is only 
weakly affected by the measurement of the feedback loop. For 
only a few atoms the dependence on the chosen measurement 
resolution indicates an increasing amount of atom-atom cor- 
relations with decreasing measurement uncertainty. For larger 
measurement uncertainties a a/N behavior will appear in the 
single-atom momentum variance, as would be expected classi- 
cally when the measurement uncertainty is equally distributed 
among all atoms [cf. Eq. (p"oj)]. 
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AQAV/(A Po Aq ) 




1 3 5 7 9 11 13 15 
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FIG. 3: Uncertainty product of single-atom position and momentum, 
scaled in units of the ground-state uncertainty Aq Ap , versus atom 
number N (same parameters as in Fig. 2). The minimum uncertainty 
product corresponds to the value 1. 



In Fig. [3] the uncertainty product of single-atom momentum 
and position, including now also the measurement-induced 
back-action noise in the position, is plotted versus the atom 
number for different measurement uncertainties. For small 
atom numbers atom-atom correlations can be build up, as in 
the case a/ Ap Q = 1 (filled boxes), leading to minimum val- 
ues at certain atom numbers. For larger atom numbers the 
measurement-induced noise dominates and the scaled value 
[l + (Ap /cr) 2 ] 2 is reached. The presented features are those 
of a condensate state, for other highly correlated quantum 



states one may expect even more dramatic deviations from the 
classically expected behavior. 



V. SUMMARY AND CONCLUSIONS 

In conclusion a method has been presented for treating 
the dynamics of feedback of macroscopic observables of 
a bosonic multi-mode many-atom system. The introduced 
unique correlation function obeys closed equations of mo- 
tion for the feedback operation and the free-evolution dynam- 
ics. Thus the dynamics intermitted by feedback operations at 
predefined times can be exactly described and calculated on 
a quasi single-atom level. The description contains the full 
many-atom correlations, that are generated by the feedback 
loop. In this way the dynamics of expectation values of any 
extensive property of the many-atom system, as for example 
the energy of a dilute atomic gas, can be obtained. The advan- 
tage of the method is its exactness combined with the treat- 
ment on a quasi single-atom level, vastly reducing the number 
of degrees of freedom, that would have to be employed for a 
solution of the full many-atom problem. Applications of the 
method can be found in the context of control and stochas- 
tic cooling of dilute bosonic gases Jl8[ | or feedback of a few 
atoms in an optical cavity [|l9f p0|]. 
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